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I . INTRODUCTION 


This report summarizes the research efforts of the principal investigator 
to initiate a theoretical study of silicon superlattices and their applicability 
to improved semiconductor devices. The work described here has proven to be 
very fruitful and is currently being continued by the author at this new institu 
tion, the University of Cincinnati, Additional support for this research, in 
the form of a separate NASA grant, is now pending. 

The motivation for undertaking this investigation was two-fold. First, 
there has been a great deal of recent interest in the device-application poten- 
tial of III-V semiconductor superlattices fabricated in the form of ultrathin 

o 

(-100-500A) periodically-layered heterostructures. This interest first began 
with the suggestion of Esaki and Tsu^ that the atomic-like dimensions of such 
structures could lead to interesting and useful quantum size effects and would 
create the possibility of altering familiar transport and optical properties 
over broad physical ranges. The sophisticated technique of molecular-beam- 
epitaxy has since been used to fabricate superlattices out of compatible semi- 
conducting compounds and their alloys, most notably GaAs and Ga^^^^AJl^As, where 
the lattice constants of the two constituents can be closely matched. Very re- 

cert research on the GaAs-Ga^ _^^A jIj^As system has been both intense and sue- 
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cessful. Applications to improved semiconductor lasers and electron devices 
have been achieved and many more are expected. 

The second motivating ingredient for the present work was the natural 
question of whether or not such ideas had practical relevance to silicon, still 
the most technologically important electronic material. In particular, we were 
attracted by the prospect of combining the benefits of cost-effective materials 
processing, which is offered by current silicon technology, with the improved 
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electron-device performance promised by semiconductor superlattices. The 
modified electronic structure of the semiconductor superlattice brings with It 
the possibility of controlling basic electronic properties such as electron and 
hole mobilities. Of special Interest Is the possible application of silicon 
superlattices to faster high-speed computing devices, where one factor limiting 
the speed of the device Is the relatively low carrier mobility In silicon. 

Our specific objective duHng the first year was to Investigate the 
modified electronic band structure of hypothetical silicon superlattice struc- 
tures Including, but not confined to, the heterostructure system of most direct 
Interest: SI-SI^^^^Ge^^. Toward this end, we have now studied rather extensively 
three general prototype models of Increasing complexity: 

(1} a one-dimensional superlattice In an Ideal one-dimensional model 
host: superlattice (1) 

(11) a two-dimensional superlattice In an Ideal two-dimensional host: 
superlattice (2) 

(111) a one-dimensional superlattice In a real three-dimensional silicon 
host: superlattice (3) 

In each case, we have examined In detail how the fundamental band gap, effective 
masses and related features of the electronic structure are altered from the 
host. These studies are discussed at length 1n Sec. II below. Our conclusions 
and suggestions for future directions are given in Sec. III. 


II. COMPLETED RESEARCH 


A. Electronic Structure of Model One-and Two-Dimensional Superlattices . 

The model host in superlattices (1) and (2) was chosen to simulate a 
semiconductor like silicon yet retain a large degree of physical and mathemat- 
ical simplicity. Conceptually* thi* superlattice is formed by substitutionally 
replacing host atoms by "impurity" atoms in a periodic fashion, where the 
underlying linear or square host crystal structure has a lattice constant 
a and two atoms per primative cell. The latter are separated by a distance 
d (the bond length). It is further assumed that the electron potential of the 
impurity atoms is characterized by a constant energy shift V from the host. 
Dtherwise, the nature of the impurity atoms is arbitrary and need not be con- 
sidered further. If the second atom in the primative cell is located at the 
point of highest symmetry, the geometry of superlattices (1) and (2) is speci- 
fied by only two numbers; n^, the total number of host primative cells (each 
with either 2 host or 2 impurity atoms) contained in a given (IDD) plane of 
the superlattice primative cell; and n^j, the number of those host primative 
cells which contain only impurity atoms. The geometry and electron potential 
of superlattice (1) are schematically illustrated in Fig. 1 for the case 0^*3 
and njj=l. Here 

S'n^a (1) 

is the lattice constant of the superlattice and 

Hp-atij (2) 

is the number of atoms in the superlattice primative cell. In reciprocal space, 
the size of the Brillouin zone is correspondingly shrunk by a factor of n^ and 
the boundaries of the "mini zone" for superlattice (1) are thus given by 

“1 ^ 1 * ( 3 ) 

s - s 






5 


In order to treat the electronic structure of superlattices (1) and (2), 
an appropriate LCAO (linear combination of atomic orbitals) app»^oach was 
developed. This technique was adopted both for Its basic simplicity and 
because of the Important re' jivance of the LCAO method to real tetrahedral ly> 
bonded solids like slllconj^ as discussed In Sec. II B below. In our LCAO 
description, a minimal basis set of two atomic-llke orbitals, |a>and |3>i 
was chosen. These orbitals were assumed to be orthonormal and both of s, 
symmetry. If the bond length Is a/2 In superlattice (1) andV? a/2 In super- 
lattice (2), each host atom experiences equivalent nearest-neighbor Interactions. 
Neglecting distant-neighbor couplings, this leads to five basic LCAO parameters 
for the host atoms: two diagonal energies 


Set “ |Hl«> 

(4) 

and 


eg » <b1h1$ > ; 

(5) 

and three Interatomic matrix elements 


= «xlHla'>, 

(6) 

Vgg. » <3|H|3'>, 

(7) 

and 


Vote* * <a|H|3' >, 

(3) 

where the prime denotes that the orbital is centered on 

a different site. 

the Impurity atoms, and e« are shifted to 

cx p 


4 * * h 

(9) 


and 


( 10 ) 
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while V^, , and remain unchanged. The host electronic structure 
(Vg*0) is established by a (4 x 4) secular determinant which may be solved 
analytically for the energy bands E(k). Our approach was to arbitrarily fix 
the five LCAO parameters at convenient values and study the change in the 
electronic structure as the superlattice parameters V^, n^ and nj^ were varied. 
The LCAO parameters adopted were 

£< 3 ,- 27 . ( 11 ) 

Eg - 6V. (12) 

and 


where 


''«■ ■ ''se' ' ''a«- * 


V ■ rT>^ h\^ 


03) 


(14) 


2m la^ 

is the natural unit of energy and r is a dimensionless scaling constant. These 
parameters lead to symmetrical valence and conduction bands with a direct band 
gap of magnitude E® located at the zone center (k*0) in boih superlattices (1) 
and (2). The matrix element insures coupling between the valence and 
conduction bands and a finite E^ which scales uirectly with r. The correspond- 
ing valence-and conduction-band-edge effective masses are equal in this model, 


mu * m» “ m » 


(15) 


and scale as 1/r. The host values of E® and mVm for the cases of interest 
are summarized in Table I. 

The superlattice band structure here is determined by a (2np x 2np) secular 
determinant which can only be solved numerically in general. A full computer program 
for this purpose has been developed and was used to obtain the results discussed 
below. Within the confines of the above model, the electronic structure of 
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Table I. Parameters of the model hosts In superlattices (1) and (2) for the 
choice of parameters ■ 2V, » 6V and V^, » V^g, * -V, where V « rlf 


•2 


I 


Superlattice 


e9 

m*/m 

0) 

0 

0 

0.4051 

0) 

-V 

1.66 V 

1.3841 

(2) 

-V 

0.94 V 

I.919I 
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superlattices (1) and (2) was studied as a function of the parameters V , 
Og and njj. Calculations on superlattice (1) were performed over the para- 
meter ranges 


0 < Vg/V <4.0, 

(16) 

2 < Hg < 8, 

07) 

1/8 < njj/ng <1/2. 

(18) 


The secular determinant In superlattice (1) Is (4ng x 4ng) and hence varied 
from (8 X 8) to (32 x 32) over the range of n. given In Eq. (17). Typical 
results for the superlattice energy bands are plotted 1n Fig. 2 for the case 
Og ■ 2 and n|j * 1 and five values of Vg/V. The far left-hand panel shows the 
host bands folded Into the mini zone of the superlattice. As Vg Is Increased 
from zero, the host energy bands are separated, narrowed, and pushed up In 
energy with "minigaps" formed at the mlnizone boundaries. The quantitative 
magnitude of these effects, however, is clearly band dependent. As the size 
of the minigaps Increases, the narrowed superlattice energy bands may In some 
cases being to be pushed towards one another leading to new strong Interactions. 
In Fig. 2 this can be seen to occur for the lowest lying conduction band and 
the highest lying valence band. The superlattice band gap closes continuously 
as the strength of the superlattice potential Vg is increased. As Illustrated 
in Figs. 3 and 4, the band gap closes more quickly with V when n is raised 
for fixed or when nj^/ng Is raised for fixed Og. 

The corresponding band-edge effective mass m* is also found to vary rapidly 
with Vg. Initially, m* increases sharply because both the valence and conduc- 
tion bands are flattened, but as E^ becomes smaller the bands interact strongly 
and m* decreases to a minimum value as E^-K). This behavior is shown in Fig. 5 


4 . 

A: 
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Fig. 5. 
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for the case • 2 and nj^ » 1 (and r » 1). More rapid variations occur as 
either Is Increased for fixed n|j or as n^j/Og Is Increased for fixed n^. 

The actual minimum m* achieved at band>gap closing Increases with n. and also 
slightly with nj^/Og. Thus small values of effective mass (and hence high 
carrier mobilities) can only be achieved for relatively strong superlattice 
potentials. To obtain •’>*/(m*j,Qg^ ■ 1.384) < 0.2, one requires Vg > £9*’°**»1 .66V 
In all cases studied with this model. 

Saturation of the above effects occurs rapidly as Og Increases. In the 
hypothetical limit n », one expects the band gap to first close when 


Vj • . 1 .66 V, 


OS) 


the value of the superlattice potential which would move the valence-band edge 

of the host to the position of the conduction-band edge. This value of Vg Is 

only slightly less than the one we have obtained for Og * 8, so that we expect 

little significant change In the above results for Og > 8. 

In the case of superlattice (2), we have considered structures with both 

square (SQk) and face-centered (FC) two-dimensional symmetry. The symmetry and 

dimensions of the real -space superlattices and the reciprocal -space mini zones 

are Illustrated In the left-and right-hand portions of Fig." 6, respectively. 

The effective superlattice constant Is again s * n.a. In the SRQ superlattice, 

the number of atoms per primative cell Is 

2 


"p • 2"s ■ 

while In the FC superlattice one has 

n_ * n.^. 

P s 


( 20 ) 


(21) 


The maximum fraction of impurity atoms In the primative cell is 1/4 for the 
SQR configuration, but still 1/2 for the FC configuration, as it is for super- 
lattice (1). The size of the secular determinant to be solved for the electronic 
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2 2 2 2 

structure Is now increased to (40^ x 40g ) and (Zn^ x Zn^ ), respectively. 

We have performed calculations on the SQR and PC cases of superlattice 
(2) over the parameter ranges 

0 < Vg/V < 5.0. (22) 

2 < ng < 4. (23) 

and 

1/8 < n^/ng < 1/2. (24) 

The behavior of the energy bands with varying Vg, ng and njj show basic similari- 
ties with that of superlattice (1), but a few noteworthy differences as well. 
There is again a continuous closing of the fundamental band gap as V is in- 
creased from zero in all cases considered except PC, ng » 2, n|^ - 1, where 
the gap first broadens slightly for small Vg. This behavior is shown in Pig. 

7. As in the case of superlattice (1), the rate at which approaches zero 
generally increases with ng and the fraction of impurity atoms in the primative 
cell, n|j/np. In this regard, the PC congifuration is clearly more efficient 
than the SQR one because nj^/n^ is larger for given values of ng and njj. Un- 
like superlattice (T), however, the nature of the fundamental band gap (direct 
or indirect) can change as Vg is increased. In superlattice (1), the gap is 
direct at k=0 for all values of Vg. In superlattice (2), the gap is direct at 
k=0 only for small Vg. One or more changes to either an indirect gap (indicated 
by an arrow and a numeral II in Pig. 7) or to a direct gap away from the zone 
center (indicated by an arrow and a numeral I in Fig. 7) occur in every case. 

In addition, when a band edge moves away from k=0, both longitudinal and trans- 
verse effective masses arise. Only the longitudinal masses were investigated 
in detail and these were found to depend on Vg in much the same manner as m’^ 
in superlattice (1). 
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B. Electronic Structure of One- Dimensional Silicon Superlattices 

In superlattice (3), we replaced the simple model host by a full three- 
dimensional diamond structure. The superlattice Itself was formed by replac- 
ing adjacent pairs of (100) planes of silicon atoms with Impurity atoms, so 
that the periodicity in the y and z directions is unchanged, but in the x 
direction becomes 

s - Og |, (25) 

where a is the host siliccn lattice constant which is assumed to stay fixed. 

In addition, we chose n^ to be an even integer, so that the superlattice always 
has simple tetragonal symmetry and the number of atoms per primative cell is 

Op • 2ng. (26) 

The minizone in reciprocal space is then in the shape of a square prism 
rotated 45® about the x axis, as shown in the bottom portion of Fig. 8. 

The simple minimal-basis LCAO method discussed above had to be suitable 
generalized for superlattice (3). In real covalent semiconductors like silicon, 
a minimal basis set consists of one £-like and three £-like atomic orbitals. A 
careful study of the electronic structure of pure silicon revealed that in this 
basis one must retain both nearest-neighbor and second-neighbor interatomic 
interactions to obtain a satisfactory description of the energy bands in the 
vicinity of the fundamental energy-band gap. The interatomic matrix elements 
which enter the band structure calculation can always be written as linear com- 
binations of the basic forms V^.^, and Furthermore, these 

ssa spa’ ppa PPtt 

quantities are primarily functions of bond length d =vr a/4 in the solid 
and, following Harrison, may be expressed in the form 

''Wni “ 


(27) 



Superlattice (3) 




Fig. 8. 
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where il.il' > £ or £ and the subscript m Is either a or ir. The dimensionless 
parameters were fixed by an empirical fitting of the known pure silicon 
energy bands. The best set of values that we thereby obtained are listed In 
Table II. 

As mentioned above, the bond length Is assumed fixed when the superlattice 
Is formed, so that the are unchanged. The silicon intra-atomic energies 
and Cp, on the other hand, become 

Cg + (28) 

and 




(p) 


(29) 


for the substitutional Impurity In the superlattice. Fitting of the bulk slll- 

eff« 

(s) 


con bands gaive a value £p ■ Sg “ ^*786 eV which also remains fixed. The effect 


of the superlattice potential Is thus Incorporated into two constants, V. 
and V and represents. In general, a different potential seen by s_ and 


by £ electrons. The special case 

(s) 


s s 


(p) 


(30) 


corresponds to the situation envisaged above In superlattices (1) and (2) 
where the impurity potential is rigidly raised or lowered by a constant value 
Vg. In the context of superlattice (3), we have designated this as type-I 
behavior. Another Interesting limiting case occurs for 

Vg = (31) 

in which case £ and £ electrons effectively see equal and opposite superlattice 
potentials. Detailed consideration of the electronic structure of Ge has 
suggested that this behavior, which we have designa.'-,ed as type-II, is to be 
expected in the Si-S1^_^Ge^ heterostructure superlattice. Our studies have 
further shown that a very approximate relation between the superlattice potential 
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Table II. Parameters Eq. (27) obtained from a best fit to the bulk 

silicon band gap 



’^Wm 

Nearest neighbor 

Second Neighbor 


SS0 

-1.491 

0.029 

spa 

1.308 

0.051 

ppa 

1.485 

0.351 

PPTT 

-0.276 

-0.036 



( 32 ) 


and the Ge concentration x Is 

Vj ■ 0.5x (in eV). 

Electronic structure calculations have been performed on superlattice 
(3) In both the longitudinal [100] direction and the transverse [001] direc- 
tion In type-I as well as type-II cases over the parameter ranges 


0 < Vg < 1.0 eV 

(33) 

2 < Og < 8 

(34) 


In order to maximize the effective strength of the superlattice potential In 
each case, the ratio of Impurity-layer to silicon-layer thickness was held 
constant at unity (corresponding to n^j/Og ■ 1/2), as suggested by our model 
studies. For n^ ■ 8, this represents layers of about 11A in thickness or a 

O 

superlattice period of s * 22A. Since the size of our secular determinant Is 
(8ng X 8ng), direct calculation on thicker systems would require determinants 
larger than (64 x 64). 

Typical results that we have obtained for the band structure In the 
vicinity of the fundamental energy gap are Illustrated In Fig. 9. The top 
panel shows the pure silicon bands with the characteristic Indirect gap and 
conduction band edges along the <100> directions out toward the zone bound- 
aries. The center panel shows the effect of folding the bands back Into the 
Og ■ 4 minlzone for Vg * 0. Four of the conduction-band minima remain fixed, 
but the remaining two along the [100] and [TOO] directions are displaced to 
positions near k » 0. When the suplerlattice potential Vg is Increased from 
zero, the bands split apart, the fundamental gap closes and the [100] and [TOO] 
minima are further displaced, as shown In the bottom panel of Fig. 9. For the 
case illustrated (Vg * O.SeV; type-II), the minima in fact have moved to the 
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zone center and the gap has become direct. 

As expected from our model superlattice studies, the fundamental band 
gap was found to decrease continuously as Is Increased In all cases. The 

rate of change with Is significantly faster In the type- I I case than In 
the type-I case, as shown In Fig. 10 for n^ ■ 2. This Is to be expected on 
the basis of Eqs. (28)-(31), since the ^ and £ energies are pushed together 
In the Impurity layers for the type-II superlattice. The band gap In the 
type-I case remains Indirect over the entire parameter range of Eqs. (33) and 
(34). In the type-II case, on the other hand, there Is a tendency for the gap 
to become direct (at the zone center) for large V^, as can be seen In Figs. 9- 
11. (In Figs. 10 and 11, the region between the arrows Is that for which the 
gap Is direct at k * 0). Also as In the model studies, the gap E^ was found 
to decrease more rapidly with as n^ Is Increased. This Is shown In Fig. 
n. The effect Is not as large as In superlattice (1) (compare Fig. 11 with 
Fig. 3), nor. In fact, did we find any clear evidence of saturation as n^ 
was Increased from 2 to 4 to 6 to 8. Presumably, saturation will occur, but 
evidently thicker layers must be examined to see Its onset. 

The band edge longitudinal and transverse effective masses In superlattice 
(3) display rather complicated behavior, In general, due to the Indirect nature 
of the band gap. A clearer picture emerges by examining only the effective 
masses at the zone center k ■ 0. This k point, of course, marks the valence- 
band maximum In all cases and hence is appropriate to the actual hole effective 
mass In the superlattice. But k » 0 also marks the position of the conduction- 
band minimum in type-II superlattices for large V^, as mentioned above. In 
which case this point Is directly relevant to the electron effective mass. 

At k = 0 it was found that the longitudinal mass always increases while the 
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transverse mass decreases with Increasing V^. The former behavior Is qualita- 
tively In accord with our model studies, while the latter 1s a new three-dimen- 
sional effect and one which clearly bodes well for a mobility enhancement In 
the transverse direction. Details of the variation of the longitudinal and 

it it 

transverse hole masses, m^^^ and m^^^, with are shown In Figs. 12 and 13, 

Figure 12 also Illustrates the relatively stronger behavior In the type-II case 
than the type-I case. (The arrow denoting the cusps In the type-II curve In 
Fig. 12 marks the value of beyond which the valence-band maxima moves away 
froni i; ■ 0.) As seen In Fig. 13, the Initial rate of decrease of mj^^ Increases 
with n^ and hence layer thickness. While the net effect Is not extremely large 
on the scale of the expected In the S1-S1-|_j^6ej^ superlattice, saturation Is 
again not yet In evidence and a bigger effect can be expected In thicker systems. 

(Note In Fig. 13 that the effective mass Is given in units of the bulk silicon 

* 

heavy-hole mass 



Superlattice (3) 
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III. CONCLUSIONS 

The theoretical studies described above have led to a number of useful 
findings and should greatly help focus the direction of future Investigations 
on silicon superlattices. We have found a wide range of Interesting behavior 
possible Including the desirable features of smaller fundamental energy band 
gaps and reduced effective masses. Such effects correlate strongly with both 
the chemical (I.e., the magnitude of the superlattice potential V^) and geo- 
metrical (I.e., layer thlcknes-i parameters n^ and nj^/ng) nature of the super- 
lattice. While details vary, the following general trends are In evidence: 

(1) The effective strength of superlattice varies directly with Vg, Og, 
and n^/ng. In one-dimensional silicon superlattices one can conve- 
niently choose equal layer thicknesses to maximize nj^/Og at 1/2. 
Saturation of the superlattice strength with increasing ng (for 
fixed Vg) is expected and was found to occur rapidly In the model 
superlattices studied. Ir the real silicon superlattice, however, 

O 

saturation effects were not observed for layer thicknesses up to 11 A. 

(11) The fundamental band gap of the superlattice continuously decreases 
with Increasing superlattice strength. The direct or indirect nature 
of the gap can readily change In multi -dimensioned systems. In the 
real silicon superlattice, a definite trend toward a direct gap with 
increasing superlattice strength was observed. 

(ill) The band-edge effective mass in the superlattice or longitudinal 

direction initially increases with increasing superlattice strength. 

As shown in the model superlattice studies, only for very large (and 
probably physically inaccessible) V can the trend be reversed and 
the longitudinal mass be made to decrease below its value for Vg = 0. 
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(iv) The band-edge effective mass in the transverse direction, on the 

other hand, decreases with increasing superlattice strength. In 

the real silicon superlattice studied, the effect was found to 

increase with increasing layer thickness. 

Future theoretical work on silicon superlattices will investigate thicker- 

layer silicon systems as well as variable orientations (i.e., [110] and [111] 

as well as [100]) of the superlattice. Thick-layer calculations, however, will 

require modifications in the theoretical methods used here. Promising wave- 

function-matching techniques are currently being developed in the context of 

18-20 

III-V semiconductor superlattices. Such methods should have applicability 
to silicon superlattices. 

A complementary experimental program aimed at studying Si-Si, Ge super- 

I A A 

lattices has recently been initiated at Rockwell International and preliminary 
results do show an enhanced electron mobility in the transverse direction. 

In light of our findings, it would be extremely interesting to investigate 
whether or not this arises from a decreased transverse effective mass. 

Some of work described in Sec. II B of this report was presented at the 
1981 March Meeting of the American Physical Society held in Phoenix, Arizona.^^ 
A paper summarzing this research is also currently in preparation. 


) 
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